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Abstract. We present a new mechanism of attosecond extreme-ultraviolet (XUV) 
pulses generation from a relativistic laser-driven overdense plasma surfaces in 
the wavebreaking regime. Through particle-in-cell simulations and analysis, we 
demonstrate that the observed ultrashort XUV emission for the parameters we 
considered is predominantly due to a strong plasma-density oscillation subsequent to 
wavebreaking. The coupling of the strong density variation and the transverse fields in 
the front surface layer gives rise to the transmitted emission with frequencies mainly 
around the local plasma frequency. This mechanism provides new insights into the 
scenarios of XUV generation from solid surfaces and the dynamics of laser-plasma 
interactions. 
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1. Introduction 

Laser-generated attosecond extreme-ultraviolet (XUV) pulses are of great interest as a 
powerful tool for a number of potential applications, including exploring novel ultrafast 
dynamics with unprecedented temporal resolution[T], investigating nonlinear optics in 
the XUV region[2], and as a probe for laser-plasma interactionsP, 0]. 

A promising approach to generate such XUV sources is high-order harmonics 
generation (HHG) from relativistic laser-irradiated overdense plasma surfaces[5], which 
holds the potential and advantage to obtain XUV pulses with high brilliance. Several 
radiation mechanisms have been identified, such as coherent wake emission (CWE)jni 
El El E], relativistic oscillating mirror (ROM) [101 El El El El E], and coherent 
synchrotron emission (CSE)[ini 03 E]. While CWE predominates at moderate laser 
intensities, i.e., the laser normalized vector potential o < 1, ROM and CSE are more 
efficient for highly relativistic laser intensities with a ^ 1. So far, these mechanisms 
are well understood, and quite good agreement between the theoretical predictions and 
experiments has been achieved P, [3 El El E]- 

In addition to HHG in the direction of laser reflection, transmitted XUV emission 
from the rear target surface have also been found, which is of interest for diagnostics 
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to determine the maximum plasma density. Compared to the reflected HHG pulses, it 
also has the advantage of being used directly and not subject to additional energy loss 
due to further spectral Altering. Several mechanisms have been proposed[l9], |20l EH 
E2l EHl EH ESI ESI EZl EH] , including line emission at twice of the plasma frequency due 
to inverse two-plasmon decay [H] , and emissions explained in the contexts of the CWE 
and CSE[2ZIEH|- The analysis and demonstration of the above HHG mechanisms have 
greatly enhanced our understanding of the physics of laser solids interactions and XUV 
generation processes from solids surfaces [8]. 

In this paper, we report a new mechanism of ultrafast XUV pulse generation from 
laser-irradiated plasma surfaces in the wavebreaking regime. The attosecond XUV pulse 
is generated from the front layer of the plasma and then propagates through the foil 
target, with frequencies mainly around the local plasma frequency. Through simulations 
and analysis, we identify the underlying physics is predominately due to the strong 
plasma-density oscillation in the surface layer subsequent to wavebreaking, which we call 
wavebreaking-associated transmitted emission (WTE). We also show that the emission 
is a general process for a wide range of laser and plasma parameters in the wavebreaking 
regime. Besides offering a new option to generate ultrafast XUV pulses, the radiation 
process identifled here also provides important insights into the mechanism of XUV 
generation and the dynamics of laser-plasma interactions. 

2. Simulation setup 

Both one-dimensional (ID) and 2D particle-in-cell (PIG) simulations are carried out 
using the Virtual Laser Plasma Lab (VLPL) codej2n|- We firstly present the ID results 
to investigate the radiation mechanism in detail. Here, we demonstrate the basic idea 
mainly by considering the simplest configuration of normal laser incidence and step 
plasma density profile. As such, this mechanism can be pointed out most clearly since 
it is easily distinguished from the other mechanisms. Eor example, GWE cannot play 
a role with this geometry, because it requires oblique laser incidence and short density 
gradient[H|. Besides, the emission cannot simply be attributed to ROM either, because 
ROM only occurs in the reflection direct,iou[28]. 

The incident laser is linearly polarized in 2 ;-direction, with a Gaussian temporal 
profile az{t) = exp (—t^/r^), where — cEl/{ meCOJo) is the normalized laser 
amplitude with El and lvq the laser field amplitude and the laser frequency respectively, 
e, rUe, and c are respectively the electron mass, the elementary charge and the speed 
of light in vacuum, and r = O.STq is the pulse duration with T = 2tt/u}q the laser 
pulse duration. Here for simplicity we firstly consider this quasi-single-cycle pulse. The 
effect of multi-cycle laser pulses will be discussed in section 5.2. The fully ionized 
plasma, with a thickness of d=120 nm and an electron density of Hq = lOOuc, is initially 
located between x = 5Ao and x = 5.15Ao, where Ao= 800 nm is the laser wavelength. 
He = meC(jQ/47re^ is the critical plasma density. The ions are assumed to be immobile 
due to the short interaction time being considered. The cell size is Ao/2000 and each 
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Figure 1. (a), (c) Temporal profiles and (b), (d) frequency spectra of the XUV pulses 
observed at the rear side. Frames (a) and (b) are for a/, = 20, and frames (c) and (d) 
for = 30. The insets of (a) and (c) show the temporal profiles plotted as intensity. 


cell is filled with 100 macroparticles. 

3. Radiation features 

To show the essential signatures of the transmitted XUV pulses, we present the results 
of two reference cases in Fig. [U with frames (a) and (b) for a laser amplitude = 20, 
and frames (c) and (d) for ai = 30. From the temporal profiles of the electric field 
depicted in Fig. (H^a) and (c), one can see intense few-cycle pulses have been generated. 
The emitted pulse has the same polarization with the laser pulse, i.e., only with electric 
field along ; 2 -direction. This is different from CWE, which emits y-polarized HHG even 
for a z-polarized obliquely incidence laser [8]. The peak electric field is only about one 
order of magnitude smaller than that of the laser field. The transmitted XUV emission 
is strong. Its energy is about one-fifth of that carried by the reflected XUV harmonics 
(XUV frequency components from IOcuq to IOOcuq are compared). For = 30, the 
XUV pulse is nearly single-cycle, reaching an extremely high peak field of E^IEl = 0.3, 
corresponding to E^ = 3.5 x 10^^ V/m. The insets of Fig. [11(a) and (c) show the pulse 
intensity. The full-width at half-maximum (FWHM) of the pulses are approximately 
190 attoseconds and 20 attoseconds for = 20 and = 30, respectively. 

Fig. mb) and (d) are the Fourier spectra corresponding to Fig. (c), 

respectively. The spectra display a low-frequency cutoff at the initial plasma frequency 
ojpQ = IOcuq. Thus, the pulse can be used directly as it is already filtered by the target. 
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The pulse energy for = 20 as shown in Fig. [I][b) is mostly concentrated at the 
frequencies > ojpQ. The spectrum for aL = 30 shown in Fig. HJ^d) is broader and extends 
to higher frequencies. As such, the temporal pulse width in Fig. (H^c) is much shorter 
than that in Fig. [Il[a). When we further increase the incident laser amplitude to 
above 40, the spectra show lower frequencies indicating laser light transmission, as a 
result of strong target compression by the light pressure. 


4. Radiation mechanism 


In this section we focus on revealing the underlying mechanism of the XUV emission. 
Here we select the case of = 20 for demonstration, since the basic features of 
radiation and the interaction dynamics are similar for both cases, as can be seen from 
the Supplemental movies SMI and SM2[30]. 

We start with discussing the possible origins of high-frequency emission based on 
the expression for total transverse current leading to the transverse radiation field. 
Considering a more general case of laser at oblique incidence with an angle of 9, the 
expression for the radiation source, i.e., the total transverse current j.,., can be obtained 

as[H] 




e^ne{x, t) A{x, t) 
rUeCosd 7(2:, t) 


ec tan 9 


Zni{x, t) 


1 ne{x,t) 


cos 9 'y(x, t) - 




( 1 ) 


where A is the total vector potential, Z is the ion charge number, n* is the ion density, 
and By is the unit vector along y-axis. 

The second term is always along the y-direction and only occurs for oblique 
incidence with d 7 ^ 0. This is the source term responsible for the OWE mechanism |H]. 
It can emit radiation directly from the plasma oscillation rig, without the need to couple 
with the transverse laser field. This is obviously not the case in our scheme where 9 = 0 
and the emission is in the .^-direction. In addition to requiring oblique incidence, CWE 
is only possible in the presence of a density gradient, since its mechanism is basically 
the inverse process of resonance absorption. 

Under the condition of laser field at normal incidence with 9 = 0, the total 
transverse current is recast to be 


](xt) = (x t) 


( 2 ) 


In this case, high frequencies can be introduced by temporal modulation of the effective 
plasma density Ue/q and the vector potential A^, and Doppler upshifting effect. The 
ROM mechanism is dominated by the Doppler upshifting effect, and thus can only occur 
in the reflected direction [28]. In the following, we show that our observed transmitted 
emission can be mainly attributed to strong density oscillations at wavebreaking level 
coupled to transverse electric fields in the laser-illuminated surface layer. 
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Figure 2. (a)-(b) Spatially resolved spectrum of transverse electric field Ez. At 

each spatial observation point x, Fourier transform is carried out with respect to the 
temporal waveform recorded. The time interval is from t = 0 to t = 3.06To (marked 
as ti in frame (c)) for frame (a), and from t = 3.06To to t = S.OTq for frame (b). The 
vertical dashed white lines mark the initial plasma boundaries, (c) Spatial-temporal 
distribution of the electron density ne(x,t) in units of ric. Time ts = 3.10To is also 
marked, which indicates the onset of wavebreaking. Here, =20. 


^.7. Region of emission 

We firstly demonstrate that the transmitted emission originates from the laser- 
illuminated front layer of the target. To see this, we plot the spatially resolved spectra 
in Figs. [2](a)-(b), which give the information about when and where the emission occurs. 
The procedure to obtain these spectra is as following. First, we record the temporal 
profile of the transverse electric field £'^(xo, t) at the spatial position xq over a period of 
time. Next, Fourier transformation is carried out with respect to this temporal profile 
Ez{xo^t) to obtain Ez{xo^uj). We do this for each point of x in the range between 
X = 4.95Ao and x = 5.20Ao. Finally we map the spatial-spectra distribution of Ez{x^uj) 
as Figs. [2](a)-(b). The time interval for Fig. [2](a) is from t = 0 to t = 3.06To, and for 
Fig. [2](b) from t = 3.06To to t = S.OTq. Here Tq is the laser period. Time t = 3.06To 
is also marked as ti in Fig. [2](c), which shows the spatial-temporal distribution of the 
electron density ne(x, t). Time ts = 3.10To marked in Fig.[2](c) indicates, as we will show 
later, the onset time of wavebreaking. This means that Fig. [2](a) and (b) are plotted 
respectively before and after the time when wavebreaking occurs. As can be seen, there 
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Figure 3. (a)-(f) Spatial profiles of rie (red), longitudinal electric field Ex (black), 
transverse electric field (magenta), velocity in x-direction Vx (blue), and x- 
component of the total force Ex (green) at six reference times, to = 2.66To, ti = 3.06To, 
t 2 = 3.09To, ts = 3 .IOT 0 , t 4 = 3 .I 6 T 0 = to + To/2, and to = 3.18To, where Tq is the 
laser period. Here Ul = 20. The very sharp high-density spike formed at tz and the 
onset of a multi-stream motion in the phase space are the signatures of wavebreaking. 


is no transmitted emission before the wavbreaking occurs (see Fig. [21(a)), while after the 
onset of wavebreaking the transmitted emission is observed (see Fig. 121(b)). Note that 
the frequency of the transmitted emission in Fig. 121(b) accords with the XUV spectrum 
shown in Fig. [11(b). In addition to the timing of the transmitted emission, another 
important observation is that this XUV pulse observed at the target rear side originates 
from the front (laser-illuminated) layer of the target. The XUV pulse then propagates 
through the plasma slab. 

Next, we show that strong plasma density oscillations occur in this front layer. 
From the spatial-temporal distribution of the electron density shown in Fig. [21(c), one 
can see that electrons at the front plasma surface are initially pushed forward by the 
laser light pressure and then bounce back in the first half laser cycle. The same process 
repeats in the second half laser cycle, but with some different features: in addition to 
follow the driving-laser-pulse shape, the plasma surface also exhibits higher frequency 
oscillations. As mentioned above, it is during this time the transmitted emission occurs. 

4 . 2 . Onset of wavebreaking 

To understand how the strong density oscillations arise, we present the plasma dynamics 
at six characteristic times in Fig. El to = 2.66To, ti = 3.06To, ^2 = 3.09To, ts = 3.10To, 
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t4 = 3 .I 6 T 0 = to + Tq/2, and = S.ISTq. In the first half laser cycle, due to the 
large laser ponderomotive force, surface electrons are pushed deep inside the plasma 
(see Fig. [31(a)), creating a large electrostatic field. As the laser ponderomotive pressure 
passes its first maximum, the surface electrons are pulled back by the large electrostatic 
restoring force, gaining a large kinetic energy before exiting beyond the initial foil edge 
at X = 5Ao (see Fig. 131(b)). This group of electrons experiences a stronger inward 
acceleration in the second half of the laser cycle, when the ponderomotive pressure and 
electrostatic forces in the vacuum region are co-directed. When this group of electrons 
returns to the plasma edge, it meets background electrons that were initially deeper 
inside the surface layer and are now moving in the opposite direction. When the two 
groups of electrons cross (see Fig. 131(c)), a very sharp high-density spike forms at ts 
(see Fig. 131(d)), also indicating the onset of a multi-stream motion in the phase space - 
signatures of wavebreaking |3 1 1 . For initially cold plasma, the wavebreaking causes an 
extremely high spike in the local plasma density, although thermal pressure effects may 
limit the actual density increase [321133] • 

4.3. Strong density oscillation subsequent to wavebreaking 

The onset of wavebreaking is followed by a high level of plasma density oscillation. 
This can be seen from the spatial-temporal profile of the electron density shown in 
Fig. 121(c), and more clearly from the movie SMI in the Supplemental Material [3(1]. 
The dynamic process of the density oscillation subsequent to wavebreaking can be 
understood as following. The electrons, which start to move in the negative direction 
at the beginning of the second half laser period, largely affect the following plasma 
oscillation. They represent the electron boundary after crossing over the electrons 
moving to the right. Afterwards, the laser ponderomotive force dominates and it reverses 
the boundary electrons (see Fig. [31(e)). At the same time, the other electron bunches 
move further inside the plasma. This leads to the density profile largely different from 
that of half laser period ago when no wavebreaking occurs (see Fig. [31(a)). In the case 
of without wavebreaking, decreases exponentially in the surface layer after reaching 
its maximum, while in the case of wavebreaking the electric field decays much slowly 
(see Fig. [31(e)). As a result, the electrostatic force starts to dominate in the surface 
layer (see the total force) and consequently the electrons in this area start to move in 
the negative direction. At the same time, the electrons near the boundary still move in 
the positive direction, so a second sharp density spike is expected (see Fig. [31(f))- This 
process repeats several times that several periods of density oscillation can be expected. 

From the point of view of wavebreaking, a direct consequence of wavebreaking is a 
large number of formerly nonresonant main body electrons can rapidly exchange energy 
with the plasma wave and acquire momenta efficiently. In other words, wavebreaking 
implies a considerable fraction of the plasma electrons are trapped and start oscillating. 
This leads to a greatly increased plasma fluctuation level. The strong oscillation of 
many surface electrons can be seen from the rotating structure of the electron phase- 
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Figure 4. (a) The temporal profiles and (b) frequency spectra of the effective plasma 
density oscillation rig/y, recorded at a fixed position x = 5.05Ao near the front surface. 
The blue lines correspond to a time interval between t = 2.5To and t = S.OTq (i.e., 
before the wavebreaking), while the red lines between t = S.OTq and t = 3.5To (i.e., 
after the wavebreaking). 


space distribution at latter times, as shown by the blue dotted lines in Figs. |3t^f)-(e) and 
more clearly from the supplemental movie SMl|3nj. 

The density variation is low and smooth before the wavebreaking, while high and 
fast after it. For further demonstration, we present one example of the density variation 
recorded at a fixed position near the front surface x = 5.05Ao. Figures HU^a) and (b) 
show the temporal profiles and the corresponding Fourier spectra, respectively. Here 
the (relativistically corrected) effective plasma density Ug/y is used and the Lorentz 
y-factor is cell-averaged (y =< y >ceii)- The blue lines in Fig. H] correspond to a time 
interval between t — 2.5To and t = S.OTq (i.e., before the wavebreaking), while the red 
lines between t = S.OTq and t = 3.5To (i.e., after the wavebreaking). On the one hand, 
the amplitude of the density variation is much higher after the wavebreaking than before 
it (see the temporal profiles shown in Fig. 111(a)). On the other hand, the frequency of 
the density oscillation is also higher after the wavebreaking, as shown in the frequency 
spectra of Fig. 111(b). 


4.4- Emission of XUV pulses 


4.4-1. Radiation source Although the plasma oscillation is longitudinal, it is coupled to 
the electromagnetic emission via the transverse velocity of electrons in the front surface 
layer. To understand this coupling, we consider the wave equation with source term (as 
shown in equation (|2])) in the ID case. Assuming the laser potential a(“®®^(a:,t) in the 
skin layer is large, we can obtain the expression for the small emitted wave a®: 




47re^ 

rrieC 


n. 


{x,t) nQe{x,t) 


-f(x,t) yo(a;,t) 


Aaser 


(3) 


Here, noe(x,t) and yo(x,t) are the electron density and y—factor in the skin layer 
respectively in the absence of the plasma oscillations. According to equation (Ej), the 
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emission is proportional to the laser amplitude and the amplitude of the electron plasma 
oscillation. 

The XUV amplitude increases with increasing the laser amplitude can be expected. 
To see the effect of density variation, we also did the same simulations by use of a lower 
intensity laser (o^ = 1), or using a high intensity {ul = 30) but circularly polarized laser. 
In both cases, no wavebreaking and subsequent strong density oscillation can occur. As 
a result, no transmitted emission has been observed. Only in the wavebreaking regime, 
the electron density oscillation can reach such a high level that allows this kind of 
emission efficiently generated. It is also seen from the radiation source term in equation 
(l3l) that the temporal variation of Ue/q is mainly responsible for the high-frequency 
XUV emission, since the laser field changes on a much longer time scale than the 
plasma oscillation assuming Ue ric- These demonstrate the strong density oscillation 
subsequent to wavebreaking indeed plays a dominant role in the emission process. This 
is why we call this radiation wavebreaking-associated transmitted emission. 

4-4-2- Simulated transverse eurrent The density variation itself, though playing a 
dominate role in the emission process, does not fully determine the radiation properties. 
It is the retarded transverse current distribution that contains full information about 
the radiation. The transverse electric field measured from the rear side of the target is 
given by the general equation: 

E^{t, x) ^ [ jz(t- - - —,x'\ dx' + Eiaser{t, x). (4) 

C J—oo y C J 

The first term in the right hand side (RHS) of equation dH) is the integral of the retarded 
transverse current distribution and the second term is the z-polarized laser field 
Eiaser{t,x) that propagates through vacuum. To calculate the radiation at the end of 
the simulation box, we set x = lOA and use dimensionless values. Thus we can drop the 
^-dependence of the fields and obtain: 

fXc 

Ez{t) = -27V / A it-Xc + x', x') dx' + Eiaser{t), (5) 

Jo 

with Xc = 10. Since the values of A af each point of time and space can be obtained 
from the PIC simulation results, we can calculate the integral numerically as a function 
of time. The result is shown in Fig. O 

During the first half laser cycle when no wavebreaking-associated plasma oscillation 
occurs, the value of the integral matches -Eiaser{t), i-e., the transmitted field is zero. 
During the second half laser cycle, the shape of the integral has some additional 
oscillations. Using equation (El) we get the calculated transmitted field that matches 
the result obtained directly from the PIC simulations as expected (see Fig. [Sl(b)). 

The emission can also be seen directly from the evolution of E^ obtained from the 
PIC simulations, as shown in Figs. [31(a)-(f)- E^ is an evanescent wave in the skin layer 
before the wavebreaking (see, e.g.. Fig. [31(a)). The distribution of E^ changes drastically 
after the onset time of wavebreaking, e.g., as depicted in Fig. [Hl(f), Ez arises beyond the 
skin depth and propagates through the plasma foil. From the animated version of the 
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Figure 5. (a) The numerically calculated integral given by the first term in the RHS 
of equation © (red) and the second term of the laser field Eiaser{t) (green), (b) The 
sum of the two colored plots given in frame (a) (red) and the transmitted field obtained 
directly from the PIC simulation results (green). 



x/x xlX 



8.1 8.15 8.2 8.25 8.3 

CVX 


Figure 6. (a) The current distribution Jz{t^x) as a function of time and space. The 
colored oblique lines represent the paths along which Jz{t^x) has to be integrated to 
calculate the transmitted field, (b) Spatial distribution of the retarded current density 
for five different times. Each line in frame (b) corresponds to the line of the same 
color in frame (a), (c) The transmitted fields calculated numerically using equation 
([5]). Each colored point corresponds to the current of the same color in frame (b). 


field and plasma dynamics shown in movie SMI 133). one can see more clearly that the 
transmitted emission occurs right after the onset of wavebreaking. 

We can also analyze how the retarded current density oscillation evolves with time 
and how it corresponds to each value of the transmitted held. From the green line 
in Fig. ini[b), we can see that at time t = 8.128To, the left boundary of the current 
prohle reverses its moving direction and starts to move to the right. This leads to the 
transmitted held reaching a local maximum at the same time, as shown in Fig. ISl[c). 
Similar process occurs at time t = 8.158To. The left boundary of the current prohle 
reverses its moving direction again (see the blue line in Fig. |6l(b)). Accordingly, the 
transmitted held reaches another local maximum(see Fig.|ni(c)). Obviously, this behavior 
of the retarded current corresponds to the plasma-density oscillation we discussed above. 
The density oscillation leads to a fast variation of the total current and thus gives rise 
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to radiation at higher frequencies that can propagate through the plasma. 


4-5. Theoretical analysis 


To enable the transmitted emission propagate through the plasma slab, the plasma 
should locally support higher frequency oscillations which then couple to the transverse 
electric helds. Here we give a simple model to show in principle it is possible for the local 
plasma frequency to be higher than the background plasma frequency due to density 
compression effect. Let us consider the local plasma frequency at the laser-plasma 
interface. Let the plasma be overdense with Up ^ cuq. The laser is characterized by 
its normalized vector potential a(t) = eA(t)fmeC^. Electrons oscillate transversely in 
the laser field at the laser frequency with the normalized momentum pj^ = a. The 
ponderomotive force is Fp = -Vyiai ES], where the gamma factor 7 = y/l + aP. 
The ponderomotive force compresses electron density within the skin depth. The force 
balance is then 


F = —mx — 




+ 


[ ne{x')dx' = 0 , 

Jxn 


( 6 ) 


dx Jxo 

where n, and are respectively the ion and electron density, Xq indicates the position 
of the front of the electron skin layer (may be different from the initial position Xini — 0 
due to the compression), and x is a position inside the skin layer. Taking the spatial 
derivative of equation (jEl), we find the equilibrium electron density 

9^7 


ne{x) = ni + 


dx^ 


(7) 


Thus, the equilibrium electron density within the skin layer is, as expected, higher than 
the ion density due to the compression by the laser ponderomotive force. 

For simplicity, we consider laser field as quasi-static with respect to the high plasma 
frequency. Thus, together with the electron-density compression and the acquired 
relativistic 7 -factor of electrons, the local plasma frequency changes as well. It can 
be obtained from the equation of motion dp/dt = —E and the Maxwell’s equation 
dF/dt = nep/7[2H], where E is the electric field, so that the local plasma frequency is 


LJp{x) 


ne{x) _ ojpQ +d'^^/dx'^ 

\ l{x) \ l{x) 


( 8 ) 


where ujpQ = yfiq. Because the local plasma frequency has a strong spatial dispersion, 
plasma waves excited in this region break easily. The change in the plasma frequency 
has two sources: electron density compression and increase of the electron relativistic 
7 -factor. These two effects tend to compensate each other, but this compensation is 
incomplete. 

We assume an exponential decay of the laser field in the skin layer a{x) = 
age~‘^pO)0-^s) vicinity of Xg, which is a position inside the skin layer. Considering 

within the skin layer the field amplitude can be small, we can write 7 ~ 1 -|- a^/ 2 . To 
proceed further, since u)p{x) appears in the expression for a{x), let us firstly find the 
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zeroth-order approximation of cOp near the position Xg. Now the 7 -factor is constant and 
its derivative vanishes. Then equation (jHl) gives 

^pixs) ~ ( 9 ) 

yiixs) 

To get the first-order correction, we can write a(x) ~ and set x = Xg after 

calculating the second derivative of 7 . In this case we obtain 


l^p2 — 


\ 


l{Xs) 


( 10 ) 


For the next higher-order approximations, we can insert 00^2 instead of ujpi in the 
expression for a{x) and so on in the same way. Consequently we obtain the sequence 




p n+l 


\ 


UJ, 


pO 


-h 2a^cj^ 


T{Xs) 

The limit of this sequence is given by 


lim Upn = 

n—)-oo ^ 


^pO 


n = 1 , 2 ,... 


^pO 


j{xs) - 2al 


- 0 ? 

2^s 


( 11 ) 


( 12 ) 


Therefore we arrive at the approximate local plasma frequency at the position Xg under 
the assumption of small field amplitudes within the skin layer: 




1 - ^O? 
2^s 


(13) 


This expression is larger than the background plasma frequency mainly due to density 
compression. Thus, the local plasma oscillations can excite electromagnetic waves at 
frequencies above the background plasma frequency. These waves can propagate through 
the plasma slab and exit from the rear side of the target. When the laser intensity is too 
high, plasma oscillations can have a fundamental frequency blow that of the background. 
However, if the oscillations are nonlinear, their harmonics can propagate through. 


5. Discussions 

5.1. Parametric study 

Here we present a systematic study to show how the parameters such as laser intensity, 
initial plasma density and target thickness influence the generated XUV pulses. Figure 
EKa) shows the energy conversion efficiency rj grows with the normalized laser amplitude 
ul, reaching about 4 x 10“^ when = 30. Here, r] = J E‘^{^tr)d^tr/ f E^(^o)d^o, with 
^ = t — x/c, ^tr and denoting the transmitted and incident pulses, respectively. When 
further increases, rj decreases and the foil becomes more transparent to the laser pulse. 
Figure [7](b) suggests that higher laser intensity also favors the generation of XUV pulse 
with shorter FWHM duration At, and correspondingly broader FWHM bandwidth Acu. 
Figure [Z](c)-(d) show the influence of the initial plasma density rig. While At increases 
and Aa; decreases with the increase of Ug, the conversion efficiency rj drops rapidly. This 
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Figure 7. Influences of laser normalized amplitude ol (a-b), initial plasma density 
He (c-d), and foil thickness d (e-f) on the XUV energy conversion efficiency XUV 
pulse FWHM width Ar, and XUV pulse FWHM bandwidth Acj. In these simulations, 
when one parameter is varied, the other parameters are the same with that presented 
in section 2. The laser amplitude is = 20 in (c)-(f). 


is because wavebreaking and strong density oscillations are more difficult to drive for 
higher density plasmas. The laser intensity should also increase with the plasma density 
in order to keep the laser plasma dynamics the same, as indicated by the dimensionless 
similarity parameter S = n^jai^Uc, from the similarity theory in the ultrarelativistic 
regime ^ 1[36]. For the foil thickness d, it does not affect the conversion efficiency, 
as shown in Fig. me). This is understandable because the generation process occurs 
at the front layer of the target. However, Ar increases with d while Acu decreases (see 
Fig. [Zl(f)). This can be attributed to a result of dispersion when the pulses propagate 
through the plasmas. 

5.2. Effect of laser pulse duration 

Next we examine the WTE generation by considering two of the most concerns in real 
experimental cases, i.e., the effect of using multi-cycle laser pulses and presence of a 
finite plasma density gradient, to demonstrate the robustness of this mechanism. 

Figure IHl(a)-(b) show respectively the temporal profile and the Fourier spectrum of 
the WTE using a laser pulse with duration of 10 fs. The other parameters are the same 
with the case of = 20 in the above simulations. It is seen that a train of attosecond 
XUV pulses have been generated. The frequency spectrum contains finer structures than 
that in Fig.[T](b), as a result of the interference between different pulses in the attosecond 
pulse train. The effect of interference is more evident with a longer duration pulse of 30 
fs, as shown in Figs. IHl(c)-(d). Frequency components below lOwo can be attributed to 
a lowered plasma frequency of the relativistically heated foil. Using techniques such as 
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Figure 8. Temporal profile and the corresponding Fourier spectrum of the transmitted 
XUV emission for the case of using multi-cycle laser pulses. The laser pulse durations 
are 10 fs for frames (a)-(b) and 30 fs for (c)-(d). The other parameters are the same 
with that in Fig. [Ija). 


polarization gating, an isolated attosecond XUV pulse may be obtained. Nevertheless, 
these results shows the WTE mechanism also works by use of multi-cycle laser pulses. 

5.3. Effect of plasma density gradient 

Since the WTE process relies on the strong plasma-density oscillation subsequent to 
wavebreaking, a finite density ramp in the front of the plasma surface will affect the 
threshold of wavebreaking and thus the temporal and spectral structures of the WTE. 
From the cold nonrelativistic wavebreaking field Ewb = meOPp/e, which is dependent 
on the plasma frequency and thus the plasma density, we see that a density ramp can 
lower the wavebreaking threshold. Consequently, the presence of a pre-plasma allows 
the WTE to occur more easily. This can be seen by considering a longer plasma gradient 
length. Figure|9](a) shows the temporal profile of the WTE for the case of an exponential 
plasma density gradient with scale length L = Ao/5. The emitted pulse lasts a longer 
time of several femtoseconds, compared with the case of without pre-plasma in Fig. [Il[a). 
This is due to an earlier wavebreaking and the subsequent strong density oscillation 
when the laser pulses interact with the density ramp in the front of the target, as shown 
in Figs. |9ll^c)-(d). We can see the multi-stream motion of electrons from the phase 
space distribution (x, P^), indicating the wavebreaking has occurred (see Fig. El^c)). At 
the same time, the profile of the transverse field Ez shows the waveform with higher 
frequencies (see Fig. [91(d)). Compared to the case of L = Ao/5, the duration of the 
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Figure 9. (a) Temporal profiles and (b) Fourier spectra of the transmitted XUV 
emission for an exponential plasma density gradient with scale length L = Ao/5. 
(c) Phase space distribution {x^Px) and (d) profiles of the electron density Ue and 
transverse field Ez at time t = 3.4To for the case of L = Ao/5. (e) Temporal profiles 
and (f) Fourier spectra of the transmitted XUV for the case of L = Ao/10. The other 
parameters are the same with that in Fig. [U^a). 


emitted pulse is shorter for the case of a shorter scale length of L = Aq/IO, as shown in 
Fig. [9](e). Correspondingly, the bandwidth of the spectrum in Fig. [9](f) is broader than 
that in Fig. [9](b). 

5 . 4 . Multi-dimensional effects 

The results obtained so far are based on the ID case. We also performed 2D simulations 
to check whether this mechanism works in multidimensional cases. Here we only intend 
to compare with the ID results of the case shown in Fig. [la)-(b) to validate the 
basic physics. In the simulation, a very small grid step of Aq/IOOO is used in the x 
direction in order to resolve the wavebreaking related process. Each cell is filled with 8 
macroparticles. The other laser and plasma parameters are the same with those used in 
Fig.[I](^a)-(b) in the ID simulations, except that the laser pulse has a Gaussian transverse 
profile with a focal spot size of 10 Aq. Figure mil a) shows a snapshot of the transverse 
electric field distribution at time t = STq. We see that an ultrashort pulse is generated 
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Figure 10. 2D simulation results, (a) A snapshot of the transverse electric field 
distribution at simulation time t = STq. The green dashed lines mark the initial plasma 
boundaries, (b) Temporal waveform and (c) frequency spectrum of the transmitted 
emission recorded at (x=6Ao, y=0). 

at the rear side of the target. Figures ITnifhi-fe) show the temporal profiles and the 
corresponding frequency spectra of the transmitted emission observed at the position 
(x=6Ao, y=0) at the rear target side, respectively. The signatures of both the temporal 
waveform and frequency spectrum are in good agreement with the ID simulation results. 
An animated demonstration of the laser-plasma dynamics similar with the ID case can 
be seen in movie SM3 in the Supplemental MaterialjSD]- The density spike indicating the 
onset of wavebreaking at t = 3.11To and the subsequent transmitted emission at later 
times can be clearly seen from the movie. As for the concern of transverse instabilities, 
we note that experiments of femtosecond-picosecond laser interaction with nanometer- 
micrometer thin target are routinely available nowadays. A number of experiments 
have been carried out using similar parameters with ours. For example, experiments of 
a much longer pulse laser of 500 fs duration {aL ~ 20) interacting with thin targets of 
thickness 125-200 nm have successfully demonstrated the transmitted emission due to 
the CSE mechanism |18j. These results indicate the transverse instabilities are not fatal 
with the parameters we considered here. 

6. Conclusions 

In conclusion, a new regime of attosecond XUV pulses generation from over dense plasma 
surfaces, namely, wavebreaking-associated transmitted emission (WTE), has been 
demonstrated. The emission originates from the plasma front surface and propagates 
through the target, with frequencies mainly around the local plasma frequencies. The 
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underlying physics can be attributed to the coupling of the transverse fields in the skin 
layer and the strong plasma-density oscillation subsequent to wavebreaking. Thus the 
emission is evident only in the wavebreaking regime. This novel scenario of ultrafast 
XUV pulse emission from overdense plasmas provides new insights into the dynamics 
of laser-plasma interactions and the physics of radiation process. It may also offer an 
alternative option to generate relativistically intense ultrashort XUV pulses that may 
find extensive applications. 
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